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We have developed a steady state model which describes the internal cl g 
distribution of a planar dielectric sample exposed to a uniform electron beam. The 
model includes theeffects of charge deposition and ionization of the beam, separate 
mobilities tor electros and holes oleotron-ho e reco.btnatton and 
nail* production by drifting thermal electrons. If the incident beam current is 
greater than a certain critical value (which depends on sample thickness as well as 
"ample properties), the steady state solution is non-phystcal. We interpret 
this to mean that above the critical beam current, the sample breaks down. 


INTRODUCTION 


This paper describes a simple model of a beam charging experiment. The 
motivation for the model is the need to understand low voltage breakdown such as 
that which occurs in dielectric material exposed to the radiation 
space (ref. 1). Our approach to the problem is motivated by the work of 0 Dwyer 

(ref. 2) on high voltage breakdown. 

The model configuration is shown in figure 1, and it is assumed to have reached 
a steady State. A pair of infinite grounded plates are separated by an infinite, 
homogeneous dielectric, and a spatially infinite electron beam is incident ort 
arrangement normally, through one of the plates. The beam causes ionization at one 
SrraSd deposit electrons at a second rate, S. The fact that these rates are 
constant forces the solution to the problem to have symmetry about the cen ^plane 
between the plates. All variables are either symmetric ar anti-symmetric j about this 
plane. The problem considered here is simple, but sign conventions must be hand 
carefully to avoid confusion. 

We take current to be positive when it iS directed toward the right. 
Consequently, an electron beam traveling to the right represents negativ e current 
which we denote by Jr* Since the beam is losing electrons at the rate S 
n-3s-i), tl 


(electrons cm" 


the magnitude of Jg is decreasing but 

diB = eS > 0 

dx 


( 1 ) 


where e is the magnitude of the electronic charge, 
negative charge, the electric field (E) is positive 


Because of the build up of 
(i.e., directed toward the 
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right) in the left half of the dielectric (G < x < L/2), vanishes at the roidplane, 
and is negative in the right half of the dielectric. The electrons and holes 
produced by the beam drift under the influence of the electric field. We denote the 
resulting conduction current by J(> The total current is the sum of the beam 
current and the conduction current. In s today state, conservation of charge 
requires 


d... (Jr + Jc) n 0 (2) 

dx 


Since Je has the same symmetry properties as E, the solution of equations (1) 
and (2) is 


J C - (1 - 2x/L) eSL/2 


(3) 


For the geometry as we have defined it, this solution is independent of the 
ionization rate I, the mobilities of the charged species, and all other parameters 
except sample thickness and electron deposition rate. 


ELECTRON AND HOLE BEHAVIOR 


The conduction current is the sum of the electron and hole currents which are 
defined in the usual manner: 


j n = - nev n = neu n E 


XAJl 


jp = peVp = pejipE 


(5) 


where n and p are the electron and hole densities, v n and 
electron and hole drift velocities, and p n and p p are the 
mobilities of electrons and holes. The electron and hole 
symmetry properties as J(^ and E. Because of the symmetry 
will consider only the left half of the dielectric (0 < x 
only Jg, v n , and p (the net charge density) are negative, 
jp, n, and p are all positive. 


Vp are the 
trap modulated 
currents have the same 
of the problem, we 
< L/2). In this case, 
E, j n » 


With the use of equations (1) - (5), one variable can be expressed in terms of 
the others. Solving for p results in 


P - - M n n/p p + (1 - 2x/L)SL/(2Ep p ) > 0 


( 6 ) 


Transport equations can be written for electrons and holes: 


'• A t 


dUhVp) 

dx 


vn + (I + S) - knp 


lL(P v n) R »n t I - knp 
dx *■ 


wftme k in the recombination rate and v is the collision ionisation coefficient. 
Because of equation (6) the equation for holes (eq. fi) is redundant. There is no 

oWr d0ri f? UV ? h0crtu f° of the assumption of steady state, so the gradient of the 
elect* an flux is equal to the three source and sink terms on the right. Note that 
each unsigned term on the right-hand side of equation (7) is positive: the first 
term is the avalanche (i.e., collision ionization) term, the second is the bei 
ionization and charge deposition, and the third is recombination. 


The only other equation needed is that for the electric field: 


dli/dx 


e(p - n)e 


“'T(x? l S 1 t ?^ d l e i^ r t C C wfT l f the r? 1 * “ terial - Th » s ^"y forces E(x) 
JnV x / - JpU; - u at X » L/2 which arc the boundary conditions. 


the elS?on q current: (4) ~ (?> ’ ^ ° btain ^ followin ^ e <l uatian *°r the behavior of 


-dj n /dx = +vj n /(M n E) + (i + S)e 
+ Un/( e M n E)] 2 eku n /up 
- (1 - 2x/L) kSLj n /(2 Mn M p E2) 


vari ^ les in this equation are negative. Therefore, in the region under 

them d Th ?"V a11 teri " S °” the ri 8 ht-hand side have the sign explicit in front of 
them. The last two terms in combination are negative as in equation (6) The 

equation governing the electric field is equation (9). Equation can'agS be 
used to eliminate p and equation (4) to eliminate n with the result that 8 


dE/dx = -tO +Mp/nn) Jh^ C^pE) 

“ (1 - 2x/L) eSL/(2 M pE)]/e 


We have adopted the form of the collision ionization term given by reference 3: 


V » a 0 M n |e exp(-E 0 /|l-:j) 


where a Q has units of inverse length. 


W*«ei W*-. '*V' 


(L. , 
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Equations 
tor j n and K. 
dimensionless 

(10) and (ll) a*e a pair of coupled ordinar*' 
they may be cast into dimensionless form by 
variables 

differential equations 
the use of the 
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1 - 2x/L 

(13a) 


^ * 

^j n /oSL 

(13b) 


n *» 

E/aE 0 

(13c) 


and the dimcnsiunicur.; ;•••• t.-nti. 


5 eSt°/(4eppE 0 2 ) 

(14a) 

*■ el'/epn 

(14b) 

0 'p/M n 

(14c) 

A a 0 L/2 

( 14d) 

° = I/S 

( 14e) 

The two equations (10) and (11) become 


d*/du_ = Alexp(-l/a|nl) InL/n + (1 + a) - hrfn _ f ) / n 2 

(15) 

dq/du = [(1 + 6)5 - uj/ n 

(16) 


The boundary conditions are {(()) » n(0) = 0. 


PKOFEKTIES OF THE SOLUTIONS 


Note that (u - O/n is the dimensionless hole density while j/q is the 
dimensionless electron density. The requirement that these two densities be 
non-negative places a constraint on 


Any solution which falls outside this range is not physically meaningful. 

Let us first consider the case for which there is no collision ionization 
(A = 0). Then the two equations (15) and (16) hav> the solution 


t(u) = cu 


(18) 


n(u) = gu 


where c is the solution of 


(1 + 6 - b) c 2 - 1(1 + *)(1 + o) + 1 - b] c + (1 + O) = 0 


(19) 


and g is defined by 

g 2 * 1(1 + 6) c - 1] 


( 20 ) 


If the solution is to be physically meaningful, then c must he less than unity For 
many situations of interest, both 6 and b are very small so that c - 1 - * where A__ 
is much less than unity _ To second order in the srall quantities 6 and b 


6(1 - 6 - b/o) 


(21a) 


6b/o 


(21b) 


For purposes of illustration, however, it is preferable to use values of order one, 
In figure 2 wc show the solution for d= 0.5 and b = o = 1.0. Jhe solu 
independent of a, the dimensionless electron deposition rate. Both the electric 
field and the two currents (electron and hole) are linear. The n T ± s 
densities are constant. The quantity <t> is the dimensionless potential and is 

defined as 


c£(u) = 2V(x)/(aE 0 I 


u 

L) = j n(u 

0 


)du - d> c 


( 22 ) 


whore V(x) U the eloctroatatic potential with V(0) - V(L> - 0, and is the 
dimensionless potential at u = 0 (x = L/2). 

tn the preceding example the dimensionless electron and hole currents 
and u respectively) as well as the dimensionless e ^ ?£ 

independent of the dimensionless electron deposition rate (a). ( » 


dimensional quantities are strongly dependent on the electron deposition rate.) 
However, when A is non-zerc), the collision ionization term introduces an explicit 
dependence on a. Because this term is always positive, its presence causes both 
dydu and Uu) to increase. We have solved equations (15) and (16) numerically for 
tne same values of &, b, and a as above, but with A = 0.01. As expected the 
solution is no longer independent of a. Solutions for three values of a are shown 
m figure 3. Note that for a = 4.58 both the hole current and the hole density 
vanish at the electrodes. If a is increased beyond 4.58, the solution predicts 
negative hole densities near the electrode. Since negative hole densities are 
physically meaningless, this means that there is no steady state solution for 
a > 4.58. We interpret this to mean that the dielectric will break down. 


DISCUSSION 


We have presented a simple model of the effects of an electron beam on a 
dielectric sample. We have assumed that the sample is homogeneous and Lhat the 
incident beam is spatially uniform. We have also assumed-that the beam deposits 
electrons uniformly throughout the sample. We have found that if the incident beam 
current (or electron deposition rate) becomes larger than a critical value, there 
are no steady state solutions, which we interpret to be an indication of breakdown. 


We have only begun to explore the properties of this model for realistic values 
of the model parameters. We anticipate that the simplicity of the model will limit 
the accuracy with which it represents a real dielectric charging problem. However, 
we hope that the very simplicity of the model will make it possible to thoroughly 
study and understand the physical processes leading to breakdown in this idealized 
case. We feel that this is an important first step in the development of more 
realistic models which take into account material inhomogeneities (e.g., localized 
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Figure 1. - Model Geometry. Electrons moving to the right produce negative 
current. The beam deposits electrons uniformly throughout the sample. 

These electrons drift toward the two electrodes where they return to ground. 
In steady state the total current flowing to ground is equal to the differ- 
ence between the transmitted beam current and the incident beam current. 



Fiqure 2. - Solution for nd collision ionization (x = 0), 6 = 0.5, b a a = 1.0. 
Total conduction current is proportional to u. The electron current (t) and 
the hole current (u - c) are constant fractions of the total conduction 
current. The electron and hole densities [c/n and (u - e)/n] are also 
constant, n is the dimensionless electric field, and 4 is the dimension- 
less potential. 
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Figure 3. - Solutions for x = 0.01, « = 0.5, b = o = 1.0, and a 
4.58. When the dimensionless electron deposition rate a >4.58 
gives negative hole densities [(u - c)/n] near the electrodes, 
that there is no physically meaningful steady state solution, 
dielectric breaks down. 


1* 3.5, and 
the solution 
This implies 
.e., the 
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